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ABSTRACT 

IGR J16479-4514 is a Supergiant Fast X-ray Transient (SFXT), a new class of High Mass X-ray Binaries, 
whose number is rapidly growing thanks to the INTEGRAL observations of the Galactic plane. It was regularly 
monitored with Swift/XRT since November 2007, to study the quiescent emission, the outburst properties and 
their recurrence. A new bright outburst, reaching fluxes above 10~ 9 erg cirT 2 s , was caught by the Swift/BAT. 
Swift immediately re-pointed at the target with the narrow-field instruments so that, for the first time, an outburst 
from a SFXT where a periodicity in the outburst recurrence is unknown could be observed simultaneously in the 
0.2-150keV energy band. The X-ray emission is highly variable and spans almost four orders of magnitude in 
count rate during the Swift/XRT observations covering a few days before and after the bright peak. The X-ray 
spectrum in outburst is hard and highly absorbed. The power-law fit resulted in a photon index of 0.98+0.07, 
and in an absorbing column density of ~5xl0 22 crrT 2 . These observations demonstrate that in this source 
(similarly to what was observed during the 2007 outburst from the periodic SFXT IGR Jl 1215-5952), the 
accretion phase lasts much longer than a few hours. 

Subject headings: X-rays: individual: IGR J16479-4514. Stars: individual: IGR J 16479 -45 14 



1. INTRODUCTION 

Supergiant Fast X-ray Transients (SFXT) are a new class 
of High Mass X-ray Binaries associated with blue super- 
giant companions, several members of which were discovered 
thanks to the INT EGRAL observations of the Galactic plane 
( Sgue raet al.ll2005l) . They are sources with transient X-ray 
emission concentrated in short and bright flares (with a typi- 
cal duration of a few hours), a peak luminosity in the range of 
10 36 -10 37 erg s ~' and a quiescent level of 10 32 erg s _1 (e.g. 
lin'tZandl20 05). The short duration flari ng activity is part of a 
longer accretion phase at a lower level (Ro mano et al.l l2007). 

IGR J16479-4514 is a hard X-ray transient dis covered by 
INTEGRAL in 2003 August dMolkov et al.1 120031 Hard X- 
ray activity was observed on August 8 and 9, at a level of 
~ 12mCrab (18-25keV), while the following day the flux in- 
creased by a factor of ~2. Other outbursts caugh t with INTE- 
GRAL in 2003 and 2004-2005 were reported by lSguera et"ai1 
(2005) and Sguera et al. (2006), respectively, with peak fluxes 
above 10~ 9 erg cirT 2 s _1 (20-60 keV). The recurrent and short 
outbursts observed from this source led to a suggestion that it 
belongs to the SFXT class. 

A frequent hard X-ray (E>20 keV) flaring act ivity was re- 
cently discussed by Walter & Zurita Hera s (2007), who report 
on 27 short (duration <15 ks) flares and on 11 long (>15 ks) 
flares in archival INTEGRAL data, spanning times from 2003 
January 1 1 to 2005 December 2, with variable fluxes. 

IGR J 16479 -45 1 4 was observed once with XMM-Newton 
(IWalter et al.l 120061) in 2004 March, when it displayed a 
low level X-ray emission. The joint EPIC pn and IS- 
GRI/INTEGRAL spectrum was successfully fit with an ab- 
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sorbed Comptonized spectrum, with an intrinsic column den- 
sity of (7.7±1.7)xl0 22 cm~ 2 , an electron temperature kT e > 
13keV, and optical depth r < 1.8. The unabsorbed 2-100 
keV flux was 1.8xl0~ 10 erg cirT 2 s _1 . An upper limit to the 
presence of an iron line could be placed at EW< 280 eV. 

This source is normally detected in the Swift/BAT transient 
monitor 7 (15-50keV) at a level of 4mCrab and showed 23 
flares above 300 mCrab and 4.5 cr significance during the 
Swift mission up to 2008 March, includi ng the ones that trig - 
gered the Swi ft/BAT on 2005 August 30 dKennea et al.t 2005). 

2006 May 20 dMarkwardt & Krimml2006l) . June 24. and 2007 
July 29. 

As reported in iRomano et al.l (120071) . the Swift monitoring 
of the outburst of the periodic SFXT IGR Jl 1215-5952 in 

2007 February represents the deepest and most complete set 
of X-ray observations of an SFXT in outburst. We showed 
that the accretion phase during the bright outburst lasts longer 
than previously thought: days instead of hours, with only 
the brightest phase lasting less than one day. Stimulated 
by these Swift results, we are performing the first sensitive 
X-ray monitoring campaign of the activity of four SFXTs, 
IGR J16479-4514, IGR J17391-3021, IGR J17544-2619, 
and IGR J18410-0535, with the main aim of characteriz- 
ing their activity on long timescales, both the quiescent state 
and the outbursts rec urrence, and to te st our model for the 
outburst mechanism (ISidoli et al.ll2007l ). which involves the 
presence of a denser and slower equatorial wind compo- 
nent from the supergiant companion star. The results of 
the first four months of the on-going campaign with S wift, 
can be found in a companion paper (Sidol i et al.1 120081 Pa- 
per I, see Fig. [TJ, where X-ray activity outside outbursts in 
IGR J16479-4514 is discussed. Here we report on the de- 
tailed data analysis of the 2008 M arch 19 outburst caught by 
Swift/BAT (Ba rthelmv et al.ll2008l) and follo wed at softer en- 
ergies with Swift/XRT dRomano et al.ll2008l) . 

7 http://swift.gsfc.nasa.gov/docs/swift/results/transients/weak/IGRJ16479- 
4514/ 
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Fig. 1.— Swift/XRT (0.2-10keV) light curve of IGR J 16479-45 14 in 2008, 
background-subtracted and corrected for pile-up, P SF losses, an d vignetting. 
Data before 54525 MJD are data reported upon in Sid oli et al] <2008l) . The 
downward-pointing arrows are 3-<x upper limits. [See the electronic edition 
of the Journal for a color version of this figure.] 

2. OBSERVATIONS AND DATA ANALYSIS 

Swift/BAT triggered twice on IGR J 16479 -45 14 on 
2008 March 19, firs t at 22:44:47 UT (trigger 306829, 
lBarthelmv~ et al. 2008) while Swift had been observing it as 
part of our monitoring program, and then at 22:59:59 UT 
(trigger 306830). The spacecraft immediately slewed to the 
target, so that the narrow-field instruments started observing 
it ~ 113 s after the first trigger (the slew caused a ~ 50 s gap 
in the XRT data, and the XRT points up to ~ 60 s after the 
first BAT trigger were collected as one of our pointed obser- 
vations). 

The BAT data were analysed using the standard BAT soft- 
ware within FTOOLS (Heasoft, v.6.4). Mask-tagged BAT 
light curves were created in the standard 4 energy bands, 
(Pig. [2j, and rebinned to achieve a signal-to-noise of 5. BAT 
mask-weighted spectra were extracted over three time inter- 
vals strictly simultaneous with XRT data (Sect. 3, Fig. [3) from 
the three continuous streams of BAT data in Fig. [2] Response 
matrices were generated with batdrmgen. 

The XRT data were processed with standard procedures 
(xrtpipeline vO.11.6), filtering, and screening criteria by 
using FTOOLS. We considered both WT and PC data, and 
selected event grades 0-2 and 0-12, respectively. To account 
for the background, we also extracted events within source- 
free regions. Ancillary response files were generated with 
xrtmkarf, and they account for different extraction regions, 
vignetting, and PSF corrections. We used the latest spectral 
redistribution matrices (vOlO) in CALDB. All quoted uncer- 
tainties are given at 90% confidence level for one interesting 
parameter unless otherwise stated. 

3. RESULTS 
3.1. Light curves 

Figure □ shows the SwiftfKRT 0.2-10keV light curve of 
IGR J 16479 -45 14 throughout our 2008 monitoring program, 
background-subtracted and corrected for pile-up, PSF losses, 
and vignetting. All data in one segment were generally 
grouped in one point (with the exception of the March 19 out- 
burst, which shows up as a vertical line on the adopted scale). 
The monitoring program started on 2007 October 26 with ap- 
proximately two observations per week, but when the source 
showed signs of increased activity on March 10, the observa- 
tions became almost daily. 

Figure|2]shows the detailed light curves during the brightest 
part of the March 19 outburst in several energy bands. Fig- 
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Fig. 2.— XRT and BAT light curves of the 2008 March 19 outburst in 
units of count s and count s detector -1 , respectively. The XRT points up 
to ~ 60 s after the first BAT trigger were collected as a pointed observation 
part of our monitoring program. The gaps in the BAT data are caused by BAT 
event mode time intervals being limited to less than 600 s to reduce telemetry. 
The vertical dot-dashed lines in the bottom panel mark the two BAT triggers. 
[See the electronic edition of the Journal for a color version of this figure.] 

ure|3]shows the 4-10/0.3-4keV, 50-150/15-50keV hardness 
ratios. Fitting the 4-10/0.3-4 keV hardness ratio as a function 
of time to a constant model yields a value of 1.38 + 0.03 and 
xl = 3.46 for 8 degrees of freedom (d.o.f.). A general trend 
is observed for a spectral softening as the flux increases in the 
XRT bands. 

During the 2008 March 19 outburst the peak count rate ex- 
ceeded the one recorde d on 2005 August 30 dKennea et alJ 
12001 ISidoli etai]|2008h by a factor of ~ 5. The total (0.3- 
lOkeV) XRT light curve (Fig.[3h) shows an increase in count 
rate by a factor of ~ 5 in ~ 13 minutes, and an increase by a 
factor of ~ 350 in ~ 3 hours (comparing with the first earli- 
est point, not shown in Fig. [3). A corresponding increase is 
observed in the BAT flux on timescales of minutes. 

A timing analysis was performed on WT data (Ar = 884 s), 
after having converted the event arrival times to the Solar Sys- 
tem Barycentric frame. We searched for coherent periodicity, 
but found no evidence in the range lOms-lOOs. 

3.2. Spectra 

We extracted the mean spectrum of the brightest X-ray 
emission (obs. 00306829000) and performed a fit in the 0.3- 
lOkeV band of the WT data, which were rebinned with a 
minimum of 20 counts bin -1 to allow x 2 fitting. An ab- 
sorbed power-law model yielded an absorbing column of 
N H = (5.51+°g)x 10 22 cm -2 , a photon index F = 0.98 + 0.07, 
and xl = 0.939 (563 d.o.f., see Table [TJ. The unabsorbed 
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Fig. 3. — X-ray WT mode and BAT event mode hardness ratios during the 
2008 March 19 outburst, (a): XRT 0.3-10 keV light curve, (b): Ratio of 
the 4-10 and 0.3^ keV XRT light curves. The vertical dashed lines mark 
the five intervals for XRT time-resolved spectroscopy (Table [TJ. (c): Ratio 
of the 50-150 and 15-50 keV BAT light curves. [A], [B], and [C] mark the 
three intervals for joint XRT+BAT spectroscopy (Table0. [See the electronic 
edition of the Journal for a color version of this figure.] 
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high-energy cutoff power-law model (cutoffpl in XSPEC) 
yielded N H = (4.65+°|) x 10 22 cirT 2 , F = 0.16 + 0.48, 
E c = 7+i°, xl = 0.927 (562 d.o.f.). The Ftest probability with 
respect to the absorbed power-law model is 3.154x 10~ 3 (2.95 
cr), hence, since the cutoff energy is not well constrained, we 
favour the absorbed power-law model for the XRT data alone. 
We note however that a cutoff power-law yields a good fit to 
the joint XRT+BAT data. We also note that the derived Nh is 
in excess of the one along the line of sight, 1.87 x 10 22 cm 4 . 

The X-ray spectrum of the fainter emission (PC data, un- 
absorbed 2-10keV flux of 8.95 x 10~ n erg crrT 2 s~') is sig- 
nificantly softer: fitted using Cash statistics and adopting an 
absorbed power-law model model, we obtain T = 2.57+^ 

and N H = (5.76! 2 ^) x 10 22 cirT 2 (C-stat= 471.0 for 68.93% 
of 10 4 Monte Carlo realizations with statistics < C-stat). 

Upon examination of the best XRT hardness ratio (Fig. fSJ?), 
we accumulated spectra in five time bins which maximized 
the variations in hardness ratio and yielded ~ 2600-7500 
counts per spectrum. Table Q]reports our results. 

We extracted three XRT spectra simultaneous with the three 
BAT spectra in the time intervals 120-303 s, 541-723 s, and 
793-1092 s since the BAT trigger ([A], [B], and [C] in Fig.©, 
and performed joint fits in the 0.3-10keV and 14-150keV 
energy bands, for XRT and BAT, respectively, and applied 
an energy-dependent systematic error vector to the BAT data. 
Factors were included in the fitting to allow for normaliza- 
tion uncertainties between the two instruments, which were 
constrained to be within their usual ranges (0.9-1.1). Several 
models typically used to describe the X-ray emission from 
accreting pulsars in HMXBs were adopted. The results are 
reported in Table |2] and Fig.@] We note that the fit to the BAT 
data alone obviously results in a steep power law with T ~ 3. 

All models allow a good deconvolution of the 0.5-100 
keV emission, resulting in a rather fiat continuum below 
lOkeV together with a high-energy exponential cutoff (model 
highecut in XSPEC) at -8-20 keV (depending on the time 
interval). In the Comptonization model (compTT in XSPEC, 
ITitarchukl Q 994). we adopted a spherical geometry for the 
Comptonizing plasma, since both geometries allow an equally 
good deconvolution. 




10 

Energy (keV) 

Fig. 4. — Spectroscopy of the 2008 March 19 outburst. Top: data from 
the second BAT observation (part [C] in Fig. [3) and simultaneous XRT/WT 
data fit with an absorbed power law with a high energy cutoff. Bottom: the 
residuals of the fit (in units of standard deviations). [See the electronic edition 
of the Journal for a color version of this figure.] 



TABLE 1 
Spectral fits of XRT data. 



Spectrum Mid Time Wh 


r 


Xl (d.o.f.)/ 


(MJD) (10 22 cm" 2 ) 




C-stat (%) 



Total WT 54544.9551 
Parti 54544.9511 



Part 2 
Part 3 
Part 4 
Part 5 
Total PC 



54544.9522 
54544.9555 
54544.9569 
54544.9596 
54544.9986 



5.51+ 
6.05+ 
4.92+ 
5.07+ 
5.50+ 
5.91+ 
5.76+ 







j_zj 



+0.07 

0.83+H 

1.02 
0.93 
1.12 

2.57 



iftB 



0.939 (563) 
0.947 (95) 
0.997(149) 
0.982 (249) 
0.820 (157) 
0.796 (128) 
471.0 (68.93) a 



1 Cash statistics and percentage of Monte Carlo realizations with statistic <C-stat. 



4. DISCUSSION 

We report on the first simultaneous broad band X-ray spec- 
trum of IGR J16479-4514 in the 0.3-100 keV energy range. 
The source emission is well fit with the spectral models usu- 
ally applied to the ac creting X-ray puls ars: power laws with 
a high e nergy cutoff (|White et al.| [l983) and Comptonization 
models (ITitarchukl [19 94). The resulting parameters are also 
very similar to those of this kind of X-ray binary sources. 

The time-resolved spectroscopy of the joint XRT-BAT 
spectrum during the bright flare reveals that the initial part of 
the flare (part [A]) shows the flattest continuum below lOkeV, 
together with the lowest high energy cutoff, with respect to the 
subsequent brighter part of the emission (see Table [2] parts 
[B] and [C]), in all the three models adopted. No evidence 
is found for variability in the absorbing column density dur- 
ing the evolution of the bright flar e, and also in com parison 
with the out-of-outburst emission (Sid oli et alll200 8?>. From 
optical/IR observations of the supergiant companion, Rahoui 
et al. (2008) derived an extinction Ay = 18.5 mag, which 
translates into an absorbin g column density of 3.3 xlO 22 cuT 2 
dPredehl & Schmittl 119951) which is in excess of the total 
Galactic in the source direction (see Sect. 3). Thus there is 
evidence for neutral absorbing matter local to the binary sys- 
tem. Among the different models used to fit the broad-band 
spectrum, the Comptonization model yieded an Nh compat- 
ible with the optical extinction, thus there is no strong evi- 
dence that an excess of absorption is local to the X-ray source 
(and that could be linked to the accreting clumps, if the super- 
giant wind is inhomogeneous). This could imply that the bulk 
of the neutral absorbing matter is circum-binary matter, in- 
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stead of being associated with the accreting clumps (whose 
matter is very likely totally ionized when approaching the X- 
ray source). Although the statistics allow a constant value of 
jVh (Table |2]i, a trend is found in the hardness ratio (Fig. [3J5), 
which is decreasing with increasing flux, during the outburst. 
This is probably the indication of a reduced neutral gas col- 
umn density, as the circumstellar gas is progressively ionized 
by the X-ray flare. 

Before these observations, IGR J 1 1215-595 2 was the only 
SFXT observed in depth during an outburst dRomano et alj 
120071 ISidoli et alJ 120071) . In that case, the occurrence of 
the o utburst was predic table thanks to the periodic recur- 
rence (ISidoli et alJl2006l) . For IGR J16479-4514 the outburst 
could only be cau ght thanks to our Swift monitoring campaign 
(ISidoli et alJl2008n . 

The source light curve before 5452 5 MJD (Fig . [T]) has al- 
ready been reported and discussed in Sidoh ^etall d2008l) . It 
shows a smoothly variable flux, with a dynamic range of 
more than one order of magnitude, which apparently is si- 
nusoidally modulated with a period of ~24 days (although we 
caution that only about 2 periods of this suggested periodic- 
ity have been covered). During the most frequent monitoring 
(between 54535 and 54544 MJD) preceding the outburst, a 
highly variable source intensity was observed, with rates go- 
ing up and down on short time scales of 1-2 days and span- 
ning two orders of magnitude at maximum, before reaching 
the brightest peak. In the total light curve monitored with 
Swift, IGR J16479-4514 spanned almost four orders of mag- 
nitudes in flux in a few days. Moreover, the flux is highly 
variable on different time scales, from seconds to minutes, 
days and weeks, revealing a very i ntense flaring ac tivity, both 
during the outburst and outside it (Sidoli et al. 20081). 

The rise time to the peak (at an average rate of 
-17 counts s -1 on 54544.95 MJD) is at most ~3 hrs. Af- 
ter the peak (also monitored by the BAT), the nearest XRT 
observation (~ 1.175 days later) found the source at a much 
fainter rate (0.085+0.017 counts s _1 ). This implies that the 
brightest phase of the outburst lasts one day at most, which 
is a behaviour similar to what observed in IGR J 1 1215-5952 
(Roma no et al.ll2007l) . where the bright emission lasted less 
than 1 day, and then was followed by a much weaker emis- 
sion, although highly variable, with several shorter (but 
fainter) flares. This is also consistent with the duration of 
p revious outbursts de tected by the BAT. 

iRahoui et al.l d2008l) classify the optical companion as an 
08.51 star located at ~4.9 kpc. This, together with the 
highly variable X-ray flux, confirms the classification as a 
SFXT. In these sources, a subclass of HMXBs, the accre- 
tion onto the compact object is very likely through the strong 
wind from the supergiant donor. The orbital parameters of 
IGR J16479-4514 as well as the nature of the compact object 
are still unknown, although the X-ray spectral properties are 
similar to those of accreting pulsars. Thus, the compact ob- 
ject is probably a neutron star, as in IG R Jl 1215-5952, wh ich 
hosts a pulsar with a period of ~187 s (Sw ank et alJ l2007). 

Assuming the compact object is a neutron star and the typ- 



ical parameters for an O-type supergiant (mass, M ~ 50- 
60 M Q , a radius R = 25-30 R Q , a beta-law for the su- 
pergiant wind with an exponent /3 = 1, a wind termi- 
nal velocity of 1800 -2200 km s' 1 , and a wi nd mass loss of 
2.5xlO~ 6 M yrs -1 ; IWilson & Dopital [l985|) . we can com- 
pare the expected X-ray luminosity for Bondi-Hoyle accre- 
tion with the observed X-ray emission. The fast (< 3 hrs) rise 
to the brightest emission during outburst is difficult to explain 



TABLE 2 

Spectral fits of simultaneous XRT and BAT data. 



Spectrum 






Parameters 






HIGHECUTPL 3 




r 


£ c (keV) 


Et (keV) 


xl (dof) 


Part A 
PartB 
PartC 


5 6 +L2 

J -0.96 
5.7+ 04 

6 2+°^ 
°- -0.5 


n 74+0.29 
"■' -0.29 

1 15+ - 15 


7 2 +2 - 2 
'■ -1.5 

6.9±1.0 

6 6 +i.o 

"■"-0.7 


q Q+2.9 
17 

15 3+ 2 - 7 
1J - -1.8 


0.823 (150) 
0.982 (264) 
1.046 (301) 




£0.5-10° 


i-0.5-100 C 








Part A 
PartB 
PartC 


0.94 
2.66 
2.82 


2.05 
5.46 
5.73 








CUTOFFPL a 


Nh 


r 


E c (keV) 




xl (do!) 


Part A 
Part B 
PartC 


54+ l.2 

3 - 8 -0.5 
6.5±0.6 


0.17±0.40 
1.0±0.2 
0.97±0.17 


7 4 +23 
'■^-1.5 

15.3!| 
13.5! 2 | 




0.845 (151) 
1.042 (265) 
1.108 (302) 


coMpTT b 


N H 


To (keV) 


T c (keV) 


T 


Xl (dof) 


Part A 
PartB 
PartC 


9+1.0 
J - -0.7 
2.8±0.4 

3.2±0.4 


1 57 +0 - 24 
1-J '-0.30 

1.36±0.11 
1 41+0.12 
-0.10 


7.1+ 33 
12.0+11 
10.8+ff 


7.81 2 

5 3+ 10 

5^ 
J - -2.2 


0.829 (150) 
0.962 (264) 
1.00(301) 



a A/h is the neutral hydrogen column density (Xl0 22 cm -2 ), Y the power law photon index, 
E c the cutoff energy (keV), E{ the exponential folding energy (keV).* 3 Tq is the temperature 
of the Comptonized seed photons, r e and r are the temperature and the optical depth of 
the Comptonizing electron plasma (spherical geometry). c In units of 10 37 erg s _1 derived 
assuming a distance of 4.9 kpc. 



with an enhanced accretion rate when the compact object ap- 
proaches the O-type supergiant companion, along its orbit, if 
we assume a symmetric and homogeneous wind. The same 
behaviour was obs erved in the periodic IGR J 11215-5952 
(Sidol i et alj 120071) . In this framework, the fainter X-ray 
emission (10 34 — 10 35 erg ) preceding the bright outburst can- 
not be explained with an orbit with a period shorter than 
~ 15 days, and an eccentricity lower than ~0.5. A more de- 
tailed light curve modeling cannot be performed at this time, 
since the orbital parameters are currently unknown. 
Facilities: Swift. 

We thank the Swift Team for making these observations 
possible, particularly the duty scientists and science planners. 
PR thanks INAF-IASF Milano, where part of the work was 
carried out, for their kind hospitality. This work was sup- 
ported in Italy by MIUR grant 2005-025417 and contract 
ASI/INAF 1/023/05/0, 1/088/06/0, and at PSU by NASA con- 
tract NAS5-00136. 
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